Multiple Magnon Modes and Consequences for the Bose— Einstein Condensed Phase 
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The compound BaCuSi2C>6 is a quantum magnet with antiferromagnetic dimers of S = 1/2 
moments on a quasi-2D square lattice. We have investigated its spin dynamics by inelastic neutron 
scattering experiments on single crystals with an energy resolution considerably higher than in an 
earlier study. We observe multiple magnon modes, indicating clearly the presence of magnetically 
inequivalent dimer sites. This more complex spin Hamiltonian leads to a distinct form of magnon 
Bose-Einstein condensate (BEC) phase with a spatially modulated condensate amplitude. 

PACS numbers: 75.10.Jm; 78.70.Nx; 05.30.Jp 



The investigation of field-induced quantum phase 
transitions (QPTs) in magnetic insulators continues to 
enrich our understanding of the possible quantum ground 
states of matter [3, 0, y, 0, 0]. Structurally dimerized 
quantum spin systems play a preeminent role in these 
studies. One recent example is the phase diagram of the 
compound BaCuSi206 0,0, .8, 9], in which antiferromag- 
netic dimers formed by the S — 1/2 magnetic moments 
from pairs of Cu 2+ ions are arranged on a quasi-two- 
dimensional (2D) square lattice. This material is one of 
the best candidates for the investigation of field-induced 
BEC of magnetic quasiparticles, and a 3D-2D dimen- 
sional crossover at the QPT has been reported Q. 

The ground state of weakly interacting antiferromag- 
netic dimers is a spin singlet (\S,S Z ) — 1 0, 0) ) separated 
by an energy gap A from excited triplet states (|1, 0) and 
|1, ±1)). A QPT occurs at the field H c = A/gfi B , where 
the S z = +1 component condenses. At H > H c the av- 
erage triplet density becomes finite (magnon BEC) and 
the ground state changes from a nonmagnetic singlet to 
an ordered magnetic phase. The triplet quasiparticles 
may also be considered as hard-core bosons with a ki- 
netic energy and an effective repulsion. Depending on 
the balance between these terms, a characteristic spatial 
modulation of the triplet density may arise, whereas in a 
uniform BEC this is identical for all sites. 

The classes of field-induced QPT known to date in 
dimer spin systems are summarized in Fig. [TJ For mag- 
netic interactions with weak or no frustration, the kinetic 
energy is dominant and the ordered, or BE-condensed, 
ground state is uniform at H > Hr, as in TIC11CI3 
E 0, EE E3- In SrCu 2 (B0 3 )2 d El El triplet hop- 
ping is suppressed by geometrical frustration and the re- 
pulsion causes the condensed triplets to form a superlat- 



tice with spontaneous breaking of translational symmetry 
and the appearance of magnetization plateaus (l4j |. All 
dimer sites have finite, if weak, triplet density at H > H c 
[l3j ]. Magnetization plateaus occur also for strong and 
explicit translational symmetry breaking, which leads to 
inequivalent dimer sites and multiple magnon modes, as 



in NH4CUCI3 [15j,[16[ where the separation of the singlet- 



triplet gaps exceeds the magnon band widths. 

Prior to this study the magnetic Hamiltonian used 
to describe BaCuSi20g was based on an inelastic neu- 
tron scattering (INS) investigation with coarse energy 
resolution Q and on fits to thermodynamic data @> 8]. 
Motivated by the exotic low-temperature behavior 9J], 
the lack of high-resolution data and reports of a struc- 
tural phase transition around 100 K [l?], EH, we have 
investigated the magnetic excitation spectrum by high- 
resolution INS. We find explicitly broken translational 
symmetry in the form of structurally inequivalent dimer 
sites, which produces a BEC with spatial modulation of 
the triplet density, but without magnetization plateaus. 
As a consequence our conclusions extend beyond the de- 
termination of the spin Hamiltonian to demonstrate that 
BaCuSi206 represents a fourth category [Fig. QJIV)] of 
spatially modulated BEC system, and we propose a closer 
examination of the magnetic properties around the field- 
induced quantum critical point (QCP). 

Five single crystals of BaCuSi206 were coaligned to 
provide a total mass of 1.25 g. The total mosaic spread of 
the assembly was excellent and matched the instrumen- 
tal resolution [Fig. [Ha)]. A tetragonal-to-orthorhombic 
distortion occurring below T rs 100 K was detected very 
recently by high-resolution x-ray diffraction [18]. The 
small difference between a- and 6-axes is related to a 
structural modulation with wave vector (0 0.129 0); crys- 
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FIG. 1: Physical properties at a field-induced QPT in dimer- 
based S = 1/2 spin systems, (a) triplet dispersion at zero 
field, (b) uniform magnetization M z , (c) dimer lattice with 
triplet condensate amplitude for H > H c represented by in- 
creasing grey scale. (I) 3D systems with uniform magnon 
BEC, e.g. TlCuCl 3 . (II) Shastry-Sutherland geometry with 
magnetization plateaus and strong contrast in condensate am- 
plitude, e.g. SrCu2 (B03)2 (field in panel (c) yielding M z — 
1/3). (Ill) Same features for 3D system with strong transla- 
tional symmetry-breaking and weak interdimer interactions, 
e.g. NH 4 CuCl 3 at H C 2 < H < H s2 . (IV) Quasi-2D system 
on square lattice with weakly broken translational symmetry 
showing no magnetization plateaus, e.g. schematic model for 
BaCuSi 2 6 at H c < H < H 3 . 



tals are twinned below this transition [18[, which should 
result in spatial averaging of in-plane properties. 

INS spectra at T — 1.8 K and zero magnetic field were 
collected on the cold-neutron triple-axis spectrometer 
TASP T§|, operated with fixed final momentum kf = 1.3 
A -1 (1.5 A -1 ), focusing Pyrolytic Graphite monochro- 
mator and analyzer, cold Be filter, and open horizontal 
collimation to gain intensity. The Gaussian energy res- 
olution of 0.14(0.23) meV (FWHM) is considerably im- 
proved compared to that in Ref . [6(] . Figure [^b) shows 
the spectrum at the AF zone center Q=(q h q k qi)=(0 3) 
and can be compared directly to Fig. 2(b) of Ref. Q. 
Surprisingly, the singlet-triplet gap has a clear multi- 
peak structure over an energy range of 1 meV. This 
splitting cannot be attributed to the multi-crystal sam- 
ple [Fig. EJa)], being observed consistently in different 
instrumental configurations including flat analyzer, and 
confirmed on different instruments (below). We conclude 
that the triplet excitation spectrum of BaCuSi20g shows 
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FIG. 2: (a) Rocking curves across two Bragg peaks of 
coaligned BaCuSi2 0s sample. Horizontal line indicates in- 
strumental resolution (full width) of neutron diffractometer 
MORPHEUS, (b-d) INS spectra of BaCuSi 2 6 at T = 1.8 K 
and H = T, experimental conditions as indicated. Momen- 
tum vectors in reciprocal lattice units (r.l.u.) of tetragonal 
space group IAi/acd with a = b — 10.01 A and c = 22.47 
A [13], i.e. (qn 4) is diagonal direction on square lattice. 



intrinsically more features than originally reported. 

The INS spectrum in Fig. [^b) is described adequately 
by two Gaussian peaks. However, while the excitation at 
lower energy transfer, denoted henceforth as mode I, has 
an energy width compatible with the instrumental resolu- 
tion, the peak at higher energies is considerably broader. 
The evolution of the inelastic signal in the square-lattice 
plane is presented in Fig. [^c) , and the multi-peak struc- 
ture is observed at every point. At (0.6 4) the slope 
of the resolution ellipsoid matches that of the dispersion 
relation, giving optimal resolution conditions [Fig. [2jd)] : 
this spectrum reveals that the excitation at higher en- 
ergy transfer consists of at least two transitions, which 
are denoted as modes II and III. However, peak II re- 
mains broader in energy than expected from the instru- 
mental resolution (cf. mode I) , and therefore may consist 
of more than one mode. 

The energy dispersion of modes I— III is extracted from 
the data by least-squares fits of the resolved peaks, and 
is summarized in Figs. [3fa) and (b). A cosine disper- 
sion is observed in the square-lattice plane: modes I— III 
disperse in parallel and are well described by the lowest- 
order perturbative expression 

E a (Q) = J a - J' a (cos(Trq h + nq k ) + cos(7rq, l - wq k )) 
+ J^(cos(2irq h ) +cos(27r ft )) (1) 
+ 2J[ cos(f qi)\ cos(irq h ) - cos(nq k )\, 

where a — 1,2,3 denotes modes I— III. Parameters J a 
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are intradimer interactions, with J' a and J£ respectively 
nearest- and next-nearest-neighbor in-plane interdimer 
interactions and j£ the effective interlayer interactions. 
The measured energies have no detectable (?; -dependence 
(bandwidth below 0.05 meV), confirming previous claims 
that BaCuSi2C>6 is quasi-2D. This is a consequence of 
the frustrated geometry of the interlayer coupling [inset 
Fig. [3(b)], and the precise determination of j£ is deferred 
to a future experiment. Here we take J[ — 0, consistent 
with the fact that the other parameters in Eq. ([I]) cor- 
respond to decoupled planes and no more general model 
is required. A fit to the complete data set gives domi- 
nant intradimer interactions J% — 4.27(1), J 2 = 4.72(1), 
and J 3 = 5.04(4), with J[ = 0.49(1), J' 2 = 0.52(1), 
J'{ = -0.07(1), and J% = -0.03(1), all in meV; we take 
J3 = 5(^1 + ^2) an d J3 = ^{J'i+J'2) to approximate the 
dispersion of mode III. Thus we find in contradiction to 
Ref. Q that the only significant interdimer interactions 
are those forming the square lattice, J' a ss 0.1J a . The 
individual gaps at the AF zone center are Ai = 3.15(5), 
A 2 = 3.62(5), and A 3 = 3.94(8) meV. 

The total INS intensity of modes I— III is modu- 
lated along (0 qi) by the dimer structure factor 
\f(Q)\ 2 sin{dqi) 2 [Fig. 0(d)], with 2d = 2.85(6) A the 
average projected intradimer separation and f(Q) the 
magnetic form factor of Cu 2+ . However, fitting intensi- 
ties I and II+ni separately gives 2dj — 2.52(18) A and 
2d II+ iu = 2.99(10) A [Fig. Eld)], a difference indica- 
tive of inequivalent dimer layers, as already suggested 
by the multi-mode dispersion. However, a model with 
uniform layers of dimer oriented parallel to the c-axis 
and modulated only in this direction would predict a 
weakly decreasing INS intensity with increasing \{qh 4) |. 
The (^-dependence of the measured intensity [Fig. [3(c)] 
shows further structure which is clearly inconsistent with 
equivalent dimer sites in each plane [and thus also with 
the simple representations in Fig. QJIVc) and Eq. ([Tjl], 
We are forced to conclude that the magnetic interactions 
are also modulated within the square planes, consistent 
with what is known of the low-temperature structure 
[HI], as well as between adjacent planes. 

The dependence of the excitation spectrum on mag- 
netic fields up to H = 4 T is shown for the resolution- 
focusing point in Fig. [4] These data were collected on 
the spectrometer RITA-II [2(| with experimental condi- 
tions similar to TASP. The H = T data in Fig. 1(a) 
confirm our conclusions concerning the multi-peak exci- 
tation structure [cf. Fig. [2(d)]. At finite magnetic fields a 
redistribution of spectral weight is observed. While mode 
I is clearly split, the transitions at higher energies are 
barely resolved due to peak overlap [Figs. H(b,c)]. One 
may proceed by fixing the peak widths at the zero-field 
values and redistributing the spectral weight according 
to the relative contributions expected for Zeeman-split 
triplet modes, namely 1/4, 1/2, and 1/4 respectively for 
S z = +1, 0, and —1 components; the only fitting parame- 
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FIG. 3: (a-b) Triplet dispersion in BaCuSi2 06 at zero field 
and T = 1.8 K: (a) in-plane along Q—(qh 4), (b) perpendic- 
ular to plane with Q=(0 qi) (see inset), (c-d) Momentum- 
dependence of (scaled) INS intensity, with solid lines as de- 
scribed in text: total intensity modes I— III (black), mode I 
(red), and modes 11+111 (blue) measured with kf = 1.3 A -1 
(solid circles) and 1.5 A -1 (open circles). 



ters are the centers of the Gaussian peaks. This approach 
describes the observed finite-field spectra quite satisfac- 
torily, whereas a model with field-independent modes II 
and III cannot [dashed lines in Figs.[4(b,c)]. Fits of indi- 
vidual energies, E a , and a common g-factor are summa- 
rized in Fig.ijd): E x = 4.58(1) meV, E 2 = 5.06(1) meV, 
and E 3 = 5.39(1) meV [cf. Fig. [2(d)], while g = 2.01(4), 
in good agreement with the value reported in Ref. [8j. 
That each of the modes I— III displays individual Zeeman 
splitting into components E^'°(H) demonstrates again 
the presence of inequivalent dimer sites. Different can- 
didate mechanisms for a zero-field energy splitting, such 
as exchange anisotropy or Dzyaloshinskii-Moriya inter- 
actions [lCj, [21 1, are probably present in BaCuSi20g at 
some small energy scale [22]. However, these would pro- 



duce only three modes at finite fields and are therefore 
excluded as the origin of our observations. 

The INS results presented here identify BaCuSi206 in 
a class of low-dimensional quantum magnets which has 
not yet been considered: despite explicitly broken trans- 
lational symmetry of the spin Hamiltonian there are no 
plateaus in the uniform magnetization (Fig.[lJ Q. Both 
the presence of four dimer layers along the c-axis of the 
large unit cell and the structural distortion at T « 100 
K introduce the potential for different intra- and inter- 
dimer interactions. Our results indicate that both inter- 
layer and weak in-plane modulations are present, and 
that these are due most likely to differential dimer tilts 
away from the c-axis. A minimal magnetic model could 
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FIG. 4: (a-c) INS spectra at three magnetic-field values in 
BaCuSi 2 6 for Q=(0.6 4), T = 1.8 K, and k f = 1.3 A -1 , 
(d) Zeeman splitting of triplet modes following procedure de- 
scribed in text, (e) Schematic phase diagram around QPT; 
intensity of red shading indicates qualitative degree of c-axis 
modulation of BEC order parameter. 



be expected to generate at least 4 excitations, consistent 
with the broad, unresolved appearance of mode II. 

The inequivalent dimer sites reflected in the INS spec- 
tra result in gaps A Q which differ by less than the band 
widths of the individual triplets, in contrast to NH4CUCI3 
HH [Figs. QJIIIa,IVa)], removing the possibility of mag- 
netization plateaus [Figs.[TpIIb,IVb)]. Each of the mea- 
sured modes shows Zeeman splitting in an applied field, 
which would drive the spin system towards a QCP oc- 
curring when the lowest gap Ai vanishes at H\ = H c . 
Above H c , the density of triplet quasiparticles would be 
largest for sites of type a = 1, but because of their direct 
real-space coupling, dimers of types a = 2 and 3 would 
also acquire finite, if weak, triplet densities [3, HH . The 
BEC state must therefore be characterized by a spatially 
modulated amplitude of the order parameter. As the 
field is increased through the values H a corresponding to 
the gaps A Q (a = 2 and 3), raising the uniform magne- 
tization and the BEC order parameter, the effect of in- 
terdimer interactions in a system such as BaCuSi20g is 
to cause a mixing of the low-lying triplet modes (H) 
with the ground state. This mixing precludes further 
phase transitions and causes the magnetization to display 
a crossover rather than a kink at these fields [Fig.QJIVb)] , 
in accord with the measurement of Ref. [3] • The spatial 
modulation of the BEC order parameter may in princi- 
ple be detected by nuclear magnetic resonance (NMR) 
studies of the many inequivalent Cu sites at H > H c 

& 

In BaCuSi206 our results indicate that the dimer 
modulation occurs predominantly between planes 
[Figs. [3Jc,d)], and thus that the indices a correspond 
to separate bilayers. In this case the triplet density at 



H c < H < H2 is higher on dimer layers a = 1, and 
is lower on the intervening bilayers [cf. Fig. l(IVc)]; 
indeed the structural modulation must lead in this way 
to an enhancement of the anisotropic nature of the 
BEC order parameter, as depicted in Fig. 0Je). The 
structural modulation also raises the distinct possibility 
that geometrical frustration of interlayer triplet hopping 
[Eq. fT}] is less than perfect at low temperatures. 
However, the thermodynamic measurements of a 2D 
QCP in BaCuSi2 06 as a consequence of this frustration 
[9( demonstrate that the energy scale of any unfrustrated 
hopping component is below 30 mK. 

In summary, we have investigated the magnetic excita- 
tion spectrum of BaCuSi206 in zero field and at H < H c 
by inelastic neutron scattering on single crystals. The 
very much higher energy resolution than in earlier stud- 
ies allows the determination of a spin Hamiltonian whose 
exchange interactions indicate a complex picture of the 
ground state and of the field-induced QPT observed in 
this material. Inequivalent dimer sites give a multi-mode 
excitation spectrum, and as a further consequence the 
BEC phase at H > H c is anticipated to be spatially mod- 
ulated. Although INS studies are currently limited to 17 
T < H c , high-field magnetization and NMR experiments 
could explore in greater detail the modulated BEC phase 
whose characteristic features are predicted here. 
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